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Abstract 

The macroscopic-microscopic mass formula is further improved by considering mirror nuclei 
constraint. The rms deviation with respect to 2149 measured nuclear masses is reduced to 0.441 
MeV. The shell corrections, the deformations of nuclei, the neutron and proton drip lines, and the 
shell gaps are also investigated to test the model. The rms deviation of a-decay energies of 46 
super-heavy nuclei is reduced to 0.263 MeV. The central position of the super-heavy island could 
lie around N = 176 ~ 178 and Z = 116 ~ 120 according to the shell corrections of nuclei. 
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I. INTRODUCTION 



The concept of symmetry in physics is a very powerful tool for understanding the behavior 
of Nature. The isospin symmetry discovered by Heisenberg plays an important role in nuclear 
physics. In the absence of Coulomb interactions between the protons, a perfectly charge- 
symmetric and charge-independent nuclear force would result in the binding energies of 
mirror nuclei (i.e. nuclei with the same atomic number A but with the proton number Z and 
neutron number N interchanged) being identical [1, 2]. Although the Coulomb interaction 
can result in the isospin-symmetry-breaking (ISB), the measured energy differences in the 
excited analogue states between mirror nuclei (MED) amount to tens of keV and do not 
generally exceed 100 keV, which indicates that the "nuclear part" of the binding energies in 
pairs of mirror nuclei should be close to each other, i.e. 

E B — E c ~ E' B — E c . (1) 

Where, Eb and Ec denote the total energy and the Coulomb energy of a nucleus, respec- 
tively, and E' B and E' c denote the corresponding values of the mirror nucleus. Combining the 
macroscopic- microscopic mass formula and Eq.(l), one can obtain the constraint between 
the shell corrections of the mirror nuclei, 

|A£ - A£'| « 0, (2) 

that is to say, a small value for the difference of the shell corrections of a nucleus and of its 
mirror nucleus. It is interesting to study the constraint between mirror nuclei and the ISB 
effect for improving the nuclear mass formula, especially for the calculations of neutron-rich 
nuclei and super-heavy nuclei. 

In addition, the influence of the Coulomb interaction on the single-particle levels attracted 
a lot attention in recent years. It has been shown that single-particle effects, induced by the 
electromagnetic spin-orbit interaction and the Coulomb orbital term, produce large effects 
in the MED for nuclei in the upper sd and fp shells [3]. In [4], the authors found that the 
Coulomb potential strength does not change the position of magic gaps 50, 82 and 126, 
but strongly influences the shell structure of super-heavy nuclei. These investigations show 
that it is necessary to study the influence of the Coulomb term on the isospin-symmetry- 
breaking and on the binding energies of nuclei. The aim of the present work is to improve 
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the semi-empirical mass formula through studying the mirror nuclei constraint due to the 
isospin-symmetry and the influence of the Coulomb term on the single-particle levels and 
consequently on the shell corrections of nuclei. The paper is organized as follows: In Sec. II, 
we introduce the semi-empirical nuclear mass formula and some modifications in this work. 
In Sec. Ill, some results with the proposed model are presented. Finally, conclusions and 
discussions are contained in Sec. IV. 

II. MODIFICATIONS OF THE MASS FORMULA 

In [5], we proposed an semi-empirical nuclear mass formula based on the macroscopic- 
microscopic method [6]. The total energy of a nucleus can be calculated as a sum of the 
liquid-drop energy and the Strutinsky shell correction AE, 

E(A, Z, (3) = E LD (A, Z)U(1 + b k Pl) + AE(A, Z, (3). (3) 

k>2 

The liquid drop energy of a spherical nucleus Eu)(A, Z) is described by a modified Bethe- 
Weizsacker mass formula, 

E LD (A,Z)= a v A + a s A 2/3 + E c + a sym I 2 A + a pah A~ 1/3 6 np (4) 

with isospin asymmetry / = (N — Z)/A, and the symmetry energy coefficient, 

1 



^sym C-sym 



2 - 


|7 




2+| 


I\A 



(5) 



AV3 

The isospin dependence of the pairing term is also considered (see the expression of 5 np in 
[5] for details). The terms with bk describe the contribution of nuclear deformation to the 
macroscopic energy, and the mass dependence of bk is written as, 

2j»~ ■ (*)'"'••' *' < 6 > 

The shell correction is obtained by the traditional Strutinsky procedure [7] by setting the 
order p = 6 of the Gauss- Hermite polynomials and the smoothing parameter 7 = 1.2Huuq 
with hjjQ = 41y4 -1 / 3 MeV. For the deformation of nuclei, we only consider axially-deformed 
cases. 

In this work, we make the following modifications to the mass formula: 
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• The Coulomb energy form is slightly changed, ZiZ — 1) is replaced by Z 2 , 

Ec = a c ^- 3 [l - Z-W] , (7) 

following the form in the finite range droplet model (FRDM) [6]. This modification 
can slightly improve the rms deviation with respect to 2149 measured nuclear masses 
[8] of nuclei [N and Z > 8] by about 2 ~ 3%. 

• The microscopic shell correction of a nucleus is modified as, 

AE = Cl E sh + \I\E' sh . (8) 

Where, E s ^ and E' sh denote the shell energy of a nucleus and of its mirror nucleus, 
respectively. The additionally introduced l^l-E^ term is to empirically take into ac- 
count the mirror nuclei constraint and the isospin-symmetry-breaking effect. We find 
that this term can considerably reduce the rms deviation of masses by about 10%. 
The isospin-dependence in Eq. (8) is to consider the increase of the difference between 
neutron-neutron and proton-proton pairs in neutron-rich or proton-rich nuclei. The 
|/|-Eg h term can effectively reduce the shell correction deviation \AE — AE'\ in pairs of 
mirror nuclei, which is required from the constraint in Eq.(2) and is helpful to restore 
the isospin symmetry in the mirror nuclei. If without the |/|-Esh term in Eq.(8), we 
obtain 

\AE-AE'\=c 1 \E sh -E'J. (9) 

Considering the l^l-E^ term in AE, we obtain 

\AE-AE'\ = ( Cl ~\I\)\E sh -E' sh \ 

< Cl \E sh -E' sh \. (10) 

To illustrate this point, in Fig. 1(a) we show the values of AE' — AE between mirror 
nuclei as a function of neutron number with the WS model [5]. The balls denote the 
experimental values of the nuclear energy difference (Eb — Ec) — {E' B — E' c ) between 
mirror nuclei by adopting the Coulomb energy form in the WS model. One can see 
that the experimental binding energies of pairs of mirror nuclei are indeed close to each 
other as mentioned previously when removing the Coulomb energies and the deviations 
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FIG. 1: (Color online) (a) Shell correction difference AE' — AE in pairs of mirror nuclei. The 
squares and the open circles denote the results without and with the E' sh term being taken into 
account, respectively. The balls denote the experimental values of (Eb — Eq) — (E' B — E' c ) 
between mirror nuclei with the Coulomb energy form in [5]. (b) Deviations of the calculated 
nuclear masses from the experimental data [8]. The squares and the crosses denote the re- 
sults of WS and WS*, respectively. The solid curve denote the results of an empirical formula, 
A r = -0.7 [cos (27T ^r) + cos (2ir^)] A' 1 / 3 , which will be discussed later. 

are generally smaller than 1 MeV. The squares and the circles denote the results of 
shell correction difference without and with the |/|-Es h term being taken into account, 
respectively. Here, the shell energy of a nucleus is calculated at the deformation of 
its mirror nucleus for the sake of simplicity, since the deformations of pairs of mirror 
nuclei are close to each other for most nuclei according to the calculated results with 
WS. The WS calculations show that the shell correction differences caused by the 
Coulomb potentials are larger than 3 MeV for some mirror nuclei, which obviously 
over-predict the experimental nuclear energy differences [see the balls in Fig. 1(a)]. 
The electromagnetic spin-orbit interaction, the Coulomb orbital term or the Coulomb 
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potential strength are therefore introduced by some authors [3, 4] for improving the 
traditional Coulomb potential as mentioned previously. In this work, the influence of 
the Coulomb term is effectively considered by introducing the shell energy of the mirror 
nuclei. The shell correction difference between mirror nuclei is effectively reduced by 
about 1 MeV after the | / 1 ^ h term being considered in AE. 

• The /3q deformation of nuclei is taken into account, which slightly improves the results 
of heavy nuclei. 

III. RESULTS 

With these modifications and the obtained optimal parameters of mass formula which are 
listed in Table 1 and labelled as WS*, the rms deviations of the 2149 nuclear masses is further 
reduced by 15%, to 0.441 MeV and the rms deviations of the neutron separation energies 
of 1988 nuclei is reduced to 0.332 MeV (see Table 2). Fig.l (b) shows the deviations of 
the calculated nuclear masses from the experimental data. Considering the shell constraint 
between mirror nuclei (WS*), the results are effectively improved. 

In Fig. 2 we show the calculated shell corrections AE of nuclei with this model. Consid- 
ering the shell constraint between mirror nuclei, the nuclei with the largest shell corrections 
in the super-heavy region slightly moves to N = 176 and Z = 120. The shell energies with 
WS* for nuclei around N = 16 and N = 28 become larger in absolute values, whilst those 
for nuclei around (N = 184, Z = 82) and 100 Sn become smaller compared with the WS 
calculations. In Fig.3, we show the calculated deformations of nuclei with WS*. Obviously, 
the calculated structure of the known magic nuclei is spherical in shape. For light nuclei, the 
06 deformations of nuclei are not very obvious, compared with the intermediate and heavy 
nuclei. In Fig.4, we show the deviations of the calculated nuclear masses in this work from 
the results of other models as a function of isospin asymmetry. One can see that for highly 
neutron-rich nuclei (/ > 0.3) the deviations from these different models are large, and the 
results from FRDM and WS* are relatively close to each other, while the results from WS 
are relatively close to those of HFB-17 [10]. 

In Fig.5, we show the drip lines obtained with different mass formulas. To remove the 
fluctuations due to the shell and pairing effects, we do a polynomial fitting to the calculated 
results with the FRDM, the HFB-17, and the WS* models, respectively. The leftmost and 
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TABLE I: Model parameters of the mass formula. 



parameter 


WS 


WS* 


a v (MeV) 


-15.5841 


-15.6223 


a s (MeV) 


18.2359 


18.0571 


a c (MeV) 


0.7173 


0.7194 


c sym (MeV) 


29.2876 


29.1563 


K 


1.4492 


1.3484 


a pair (MeV) 


-5.5108 


-5.4423 


91 


0.00862 


0.00895 


fj2 


-0.4730 


-0.4632 


c\ 


0.7274 


0.6297 


V (MeV) 


-47.4784 


-46.8784 


r (fm) 


1.3840 


1.3840 


a (fm) 


0.7842 


0.7842 


Ao 


26.3163 


26.3163 



TABLE II: rms a deviations between data AME2003 [8] and predictions of several models (in MeV). 
The line a(M) refers to all the 2149 measured masses, the line a(S n ) to the 1988 measured neutron 
separation energies S n . The calculated masses with FRDM are taken from [6]. The masses with 
HFB-14 and HFB-17 are taken from [9] and [10], respectively. WS*+Ay means the correction At 
for empirically considering the tetrahedral deformation is added to the binding energy of a nucleus 
with WS*. 





FRDM 


HFB-14 HFB-17 


WS 


WS* 


WS*+A T 


a(M) 
cr{S n ) 


0.656 
0.399 


0.729 0.581 
0.598 0.506 


0.516 
0.346 


0.441 
0.332 


0.417 
0.330 



rightmost curves denote the smooth proton drip line (for odd-Z nuclei) and the smooth 
neutron drip line, respectively. The crosses denote the measured nuclei. For the proton drip 
line, the three models give similar results. For the neutron drip line, the results slightly 
deviate from each other at heavy mass region. Based on the liquid-drop model, the neutron 
separation energy of an intermediate and heavy nucleus (A ^> 1) can be approximately 
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FIG. 2: (Color online) (a) Contour plot of shell corrections of nuclei with WS*. The dot-dashed 
line passes through the areas with the known heavy magic nuclei, (b) Shell corrections of nuclei 
as a function of neutron number. The crosses and the squares denote the results of WS and WS*, 
respectively. 

written as, 

S n — —a v — 2a sym I. (11) 

Where, a v (negative value) and a sym (positive value) are the coefficients of the volume energy 
and the symmetry energy of a nucleus, respectively. For the neutron drip line (S n = 0) of 
intermediate and heavy mass region, we obtain the isospin asymmetry at the drip line 

1 a v 



t-nd 



2 flsym 



(12) 



One can see that the neutron drip line directly relates to the ratio of a v to a sym . The 
difference of the neutron drip line from different models is probably due to the difference of 
the coefficients a v and a sym adopted in the models. For nuclei with A — > oo or asymmetric 
nuclear matter (a v ~ — 16 MeV and a sym « 32 MeV), we obtain the corresponding neutron 
drip line which is also shown in Fig. 5 (solid-squared line). From the figure, one can see that 
most of measured nuclei are located in the left side of the solid-squared line. In addition, we 
show the smooth /3-stability line from FRDM (open-squared curve) and WS* (solid-circled 
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FIG. 3: (Color online) (a) Contour plot of quadrupole deformation of nuclei with WS*. (b) 
/?4 and (c) of nuclei as a function of neutron number, respectively. 
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FIG. 4: (Color online) Deviations of calculated nuclear masses in this work from the results of other 
models. The calculated masses with FRDM and HFB-17 are taken from [6] and [10], respectively. 
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FIG. 5: (Color online) Smooth drip lines from different mass formulas. The crosses denote the 
measured nuclei. The solid-squared line denote the neutron drip line of nuclei A — > oo with 
I = and taking a v = —16 MeV and a sym = 32 MeV. The open-squared and the solid- 

circled curve denote the smooth /3-stability line from FRDM and WS* calculations, respectively. 

curve) calculations, respectively. At Z = 120, the corresponding neutron number of the 
nuclei along the /3-stability line is about iV = 200. 

To further test the model, we study the shell gaps. As a measure of the discontinuity in 
the two neutron separation energy at magic neutron numbers N , the shell gap [11], 

A n (iV , Z) = S 2n (N 0l Z) - S 2n (N + 2,Z), (13) 

is a sensitive quantity to test the model. In Fig. 6, we show the calculated shell gaps at 
the magic neutron numbers No = 28, 50, 82, 126 with different models. The dashed, the 
dot-dashed and the solid curve denote the results of HFB-17, FRDM and WS*, respectively. 
The squared curve denote the experimental data. The most shell gaps can be reasonable 
well described by the WS* model, except the shell gap at sub-shell closure Z = 64 which 
is over-predicted by WS* and FRDM and is under-predicted by the HFB-17 model. In 
Fig.6(b), the peak (large shell gap) at magic number Z = 28 disappears according to the 
HFB-17 calculations, and the peak at Z = 82 can not be reasonably well described from the 
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FIG. 6: (Color online) Shell gap calculated with different models. The dashed, the dot-dashed and 
the solid curve denote the results of HFB-17, FRDM and WS*, respectively. The squared curve 
denote the experimental data. 

FRDM and HFB-17 calculations in Fig. 6(d). The experimental shell gaps at magic numbers 
Z = 20, 28, 40, 50, 82 can be remarkably well described with the proposed model. 

In addition, we study the relation between the fission barrier of super-heavy nucleus and 
the corresponding shell correction of the nucleus. Neglecting the shell energy at the saddle 
point, the fission barrier of a nucleus can be approximately written as [12], 

B f « B LD - AE. (14) 

Where Bld and AE are the macroscopic fission barrier and the shell correction of a nu- 
cleus at its ground state. For super-heavy nucleus, the macroscopic fission barrier generally 
disappears and consequently the fission barrier can be roughly evaluated through the corre- 
sponding shell correction of the nucleus. In Fig. 7, we show the fission barriers of a number 
of super-heavy nuclei (solid-circled curve) [13] which are calculated with the macroscopic- 
microscopic approach, considering the deformation of system up to Ps and the triaxial de- 
formation. The open-circled curve and the crosses denote the values of —AE with WS and 
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FIG. 7: (Color online) Fission barriers of some super-heavy nuclei. The solid-circled curve denote 
the calculated fission barriers of super-heavy nuclei in Ref.[13]. The open-circled curve and the 
crosses denote the values of — AE with WS and WS*, respectively. The dashed line denotes the 
position N = 184. 

WS*, respectively. One can see that the calculated fission barriers are generally close to 
the values of — AE with WS. Both of models in which the mirror nuclei constraint is not 
taken into account, predict the two neutron magic numbers iV = 162 and N = 178 at the 
super- heavy region. When the constraint between mirror nuclei is considered (WS*), the 
results for nuclei with Z = 116 and 118 do not change too much around N = 178, while the 
results for nuclei with Z = 114 and 120 change about 1 MeV. These investigates indicate: 
1) The calculated shell corrections (in absolute value) with the proposed mass formula are 
comparable to the fission barriers of super-heavy nuclei; 2) The mirror nuclei constraint 
could influence the shell structure of nuclei with Z = 114 and 120. In Fig. 7, the dashed 
lines denotes the position N = 184. For nuclei with N > 184, the fission barriers fall rapidly 
with the increase of neutron number. In addition, the shell corrections (in absolute value) of 
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TABLE III: a-decay energies Q a and the shell corrections in 6 a-decay chains with WS* (in MeV). 
The experimental data are taken from [14, 15]. 



A Z Q a (Expt.) Q a (WS*) A£(WS*) A Z Q a (Expt.) Q a (WS*) A£(WS*) 



90/1 
Zy4 


117 
11 1 


1 n nfi/ 1 n"\ 
!U.yO(lU J 


1 1 Q9 

11. oZ 


K 77 
-0. / / 


9QQ 

zyo 


117 
11 / 


11 1 Q/S^ 

ll.lo^o J 


1 1 £9 

11. OZ 


k 7n 
-0. / u 


zyU 


1 1 c: 
110 


1 n nn/ /in"* 

lu.uy^u j 


i n qq 
lU.oo 


-O.l / 


9QO 

zoy 


i i r: 
110 


i n /i ^/n^t 
lU.40^y J 


1 n Q/i 
lU.o4 


K 9Q 
-O.Zo 


98fi 

zoo 


1 1 Q 
llo 


Q 77/1 n~i 
it. 1 t (111 J 


y.oz 


A 


98^ 

zoo 


1 1 Q 
llo 


y.oo^y j 


lU.lo 


A AR 
-4. 40 


989 
ZOZ 


111 
111 


q 1 q /1 n~i 
y.io(iu j 


Q K9 

y.oz 


A i n 

-4. 1U 


981 
Zol 


111 
111 




i n n/i 

1U.U4 


A 1 7 
-4.1 / 


z to 


1 nn 
lUy 


n £;n/i n"\ 
y.oy(ly J 


y.00 


-4.oo 


lit 


1 nn 
lUy 




n 

y. /0 


K H7 

-0.U / 


97 A 
Z / 4 


1 D7 
1U / 


8 Q^/1 (Tl 

o.yo^iu ) 


8 71 
O. / 1 


a sn 

-4.0U 


97*? 
Z / O 


1 0.7 
1U / 




8 Q8 
o.yo 


-o.uo 


zyo 


1 on 
1ZU 




1 Q 9£ 

lo.zo 


-o.oo 


9nc 
zyo 


1 9n 
1ZU 




1 9 Q1 
IZ.Ol 


C IQ 

-O.lo 


zyz 


1 1Q 

llo 




i 9 r\a 
Iz.UO 


a 1 fi 
-0.10 


zy4 


1 1 Q 

llo 


ll.Ol^O J 


io it: 
1Z.10 


a qi 
-O.ol 


ZOO 


1 1 £. 
110 




1 1 QI 

11. ol 


£ Q9 

-O.oZ 


9nn 
zyU 


1 1 £ 
110 


1 1 nn/ £^ 
ll.UU^o J 


11 19 
11. 1Z 


-o.ol 


284 


114 




10.93 


-4.50 


286 


114 


10 33(61 


10.25 


-5.09 


280 


112 




10.90 


-4.51 


282 


112 




10.27 


-4.68 


304 


120 




12.49 


-5.08 


320 


120 




9.85 


-2.31 


300 


118 




11.70 


-5.43 


316 


118 




9.27 


-2.08 


296 


116 




10.98 


-5.43 


312 


116 




8.73 


-0.14 


292 


114 




9.12 


-5.40 


298 


114 




8.07 


-5.15 


288 


112 




9.36 


-4.46 


294 


112 




8.38 


-3.77 



nuclei with Z > 120 are obviously smaller than that of 292 1 20. According to the calculated 
shell corrections of nuclei, the central position of the super-heavy island could lie around 
AT = 176 ~ 178 and Z = 116 ~ 120. 

Furthermore, we study the a-decay energies of 46 super-heavy nuclei (the experimental 
data are taken from Ref. [14], Table I of Ref.[15] and Table II of Ref.[16]). The rms deviation 
of the a-decay energies of the 46 super-heavy nuclei falls from 0.566 MeV with FRDM to 
0.263 MeV with the WS* model (the corresponding result with WS is 0.284 MeV). In Table 
III, we list the a-decay energies Q a and the shell corrections AE in 6 a-decay chains of 
super-heavy nuclei with Z = 111 [14] and Z = 120. The available experimental data can be 
reproduced reasonably well. These calculations indicate that the proposed mass formula is 
relatively reliable for description of the masses of super-heavy nuclei. 
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IV. CONCLUSION AND DISCUSSION 



In summary, the semi-empirical mass formula based on the macroscopic-microscopic 
method has been further improved by considering the constraint between mirror nuclei. 
The rms deviation with respect to 2149 measured nuclear masses is reduced to 0.441 MeV 
and the rms deviation of the neutron separation energies of 1988 nuclei falls to 0.332 MeV. 
The shell corrections, the deformations of nuclei and the neutron and proton drip lines 
have been investigated also. The predicted central position of the super-heavy island ac- 
cording to the calculated shell corrections of nuclei could lie around N = 176 ~ 178 and 
Z — 116 ~ 120, considering the mirror nuclei constraint. The shell corrections of super- 
heavy nuclei (in absolute value) are close to the corresponding fission barriers of the nuclei 
from other macroscopic-microscopic model. The shell gaps at proton magic numbers Z = 20, 
28, 40, 50, 82 can be remarkably well described with the proposed model. The rms deviation 
of the a-decay energies of 46 super-heavy nuclei is reduced from 0.566 MeV with FRDM to 
0.263 MeV with the proposed model in this work. 

In addition, we note that the deviations from the measured masses for some nuclei with 
N « 18,26,40,56,64,70,80,88 etc. are relatively large, with both WS and WS*, which may be 
caused by the triaxial deformation of nuclei or the tetrahedral symmetry in nuclei [17, 18]. 
It is found that the strongest tetrahedral-symmetry effects appear at tetrahedral-magic 
numbers 16, 20, 32, 40, 56, etc., and the tetrahedral deformation can bring over a few MeV 
of energy gain in the nucleus [18]. We empirically describe the influence of the tetrahedral 
deformation on the binding energies of nuclei by using two cosine functions together with 
the two tetrahedral-magic numbers 16 and 20, A T = —0.7 [cos (27r^) + cos (2vr|^)] A^ 1 ^ 3 
MeV. The solid curve in Fig.l (b) denotes the results of At- With the empirical function 
At, the rms deviation of 2149 nuclear masses can be further reduced by 5%, to 0.417 MeV. 
Microscopic study on the triaxial deformation of nuclei is in progress. 
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